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Abstract 



■^ At vanishing charge carrier density, graphene on a dielectric substrate exhibits spatial doping 

^^ inhomogeneities forming electron-hole puddles [Ij. Understanding and controlling the latter is of 

1-^ crucial importance for unraveling many of graphene's fundamental properties at the Dirac point, 

,^ where the linear valence and conduction bands meet. Charged impurities trapped in the dielectric 

Q substrate [TH3] as well as the local graphene curvature [1H7] have been pointed to as the pud- 



dles' origin. Here we show by scanning tunneling microscopy/spectroscopy the coexistence and 
correlation of charge puddles and topographic ripples in graphene decoupled from a metallic Ir 



1-^ substrate. The perfect screening of charged impurities by the metallic substrate reveals charge 

Q puddles directly induced by the local graphene-substrate distance and the ensuing charge transfer. 



Furthermore, the analysis of Dirac fermion interferences shows that electron-electron interactions 
are also strongly screened. Our findings therefore demonstrate that graphene on Ir provides a new 



> 

QQ platform to study the non-interacting limit for massless Dirac electrons in graphene. 



o 



The study of factors liable to cause electron-hole puddles in graphene has so far relied on 
graphene sheets isolated by mechanical exfoliation of graphite on dielectric substrates such as 
Si02. Conflicting conclusions were drawn on the basis of different studies, some identifying 
exclusively charged impurities between graphene and Si02 as sources of the puddles [2], [3] , 
others invoking in addition the mixing of the tt and a orbitals due to local curvature [1]. 
Recent theoretical works suggest that puddles might be related in a non-trivial way to the 
rippling of graphene [5j. In this context, the strong disorder observed in graphene on Si02 
[S] and the limited knowledge about the graphene/Si02 interface calls for the use of other 
substrates. More generally, experiments based on different dielectric environments, that 
is, different strengths of charged impurities' screening have been performed. They however 
showed no significant influence of the dielectric constant on the electronic properties, thereby 
questioning the role of charged impurities for the puddle formation jOHH]. 

Metal supports are expected to fully screen charged impurities and thus are ideally suited 
for elucidating the role of other mechanisms in the puddles formation. Periodic ripples, 
arising from the lattice parameter mismatch between graphene and most transition metal 
surfaces, were for instance related to the puddles [T2] in graphene on Ru(OOOl). It however 
turned out that charge carriers in this system do not exhibit Dirac fermion-like properties, 
due to a strong hybridization between the Ad Ru and pz C orbitals [131 [II]- Even in less 
strongly coupled systems, interaction between graphene's conduction/valence bands with 
surface states of the metal [ISl US] cannot be excluded. A linear dispersion relation in 
the electronic band structure at the Brillouin zone corners was recovered in graphene on 
Ru(OOOl) intercalated with an atomic layer of oxygen below the graphene [17]. This layer 
unfortunately suppresses the graphene ripples, which prevents from addressing the possible 
relationship between puddles and topography. 

We have performed low temperature scanning tunneling microscopy and spectroscopy 
(STM/STS) on graphene prepared by chemical vapor deposition on epitaxial Ir(lll) un- 
der ultra-high vacuum (UHV). Exposure to air can decouple the graphene from its metallic 
support, transforming the moire pattern into a disordered topographic landscape (see Sup- 
plementary Information). This system thereby offers a unique configuration to study the 
impact of ripples on the electronic properties of screened graphene. 

We begin with a close look at the topography of graphene on a plain Ir atomic terrace. 
We observe a disordered topographic landscape with a typical rms roughness cr^ ?^ 100 pm 



(Fig. [Tp,), which is comparable to the corrugation of graphene on Si02. While the usual 
moire pattern is absent, atomic resolution is routinely achieved (Fig. [Ik), demonstrating 
the quality of the surface. 

The electronic density of states can then be accessed by the measurement of the local 
tunneling conductance G{V) = dl/dV. The tunnel spectra display a V-shape (Fig. ^), 
characteristic of the density of states of graphene and similar to data typically observed by 
STS on exfoliated graphene on Si02 [201 HI]. The voltage at the conductance minimum 
has been identified as the Dirac point (also called the charge neutrality point) at energy 
Ed = ^Vmin in previous works |3]. The value of Ed reflects the strength of the local 
graphene doping [Sj HI [20], which is p-type here. A dip of width of about 20 meV is also 
present around zero bias, which is frequently reported in STS studies on graphene [201422] 
and will not be discussed in this work. 

We have measured a series of high-resolution conductance maps at low temperature. At 
each position r of a tunnel conductance map, the Dirac point energy Ed is extracted using 
a parabolic fit around the minimum of G. The spatially averaged doping E^ = 340 meV 
is larger than inferred from ARPES studies in UHV, which reported Ed at 100 meV for 
graphene on Ir(lll) [15j. Maps such as shown in Fig. |2k picture the spatial inhomogeneities 
of Ed around its mean value, forming a smooth landscape of charge puddles, of about 9 nm 
size, see Supp. Info, for details. From E]^ = {hvpYT^n [2], where h and vp are the reduced 
Planck constant and the Fermi velocity respectively, we extract the distribution of the charge 
carrier density n (Fig. ^), assuming the theoretically expected value of f ^ = 1.0 x 10® m/s. 
Notably, the standard deviation cr„ ~ 9 x 10^^ cm~^ is of the same order than that reported 
in graphene exfoliated on Si02 (~ 4 x 10^^ cm~^) [3l [23], while significantly larger than 
that of graphene on a boron nitride substrate (~ 2.5 x 10^ cm~^) [23]. In the following, 
we determine the value of vp to be 0.89 x lO^m/s in our system, which leads to refining 
a„ = 1.1 X 10^2 ^^-2 

We now focus on the comparison between the Dirac point distribution and the topography. 
Since only topographic variations at length scales similar or larger than the typical puddle 
size can correlate with the charge inhomogeneities, we filter out the topographic maps from 
structures of dimensions below half the mean puddle size. Figure |2^ shows the superposition 
of an EdIy) map (color scale) along with the long wavelength-pass filtered topography z{r) 
recorded at the same position (3D profile). A very high degree of correlation between doping 



and topography is readily seen. We have quantified this by calculating the normalized cross- 
correlation function x(r) = ^^ -E_D(ri — r) x z{Y)/{aED ^ ^z) of the two data sets. The local 
cross-correlations x(0) between Ejj{v) and z{v) are in excess of 60 % in large area maps 
(Fig. pp). These correlations are independent on the region chosen, but are enhanced in 
maps with dimensions much larger than the typical puddle size. 

When correlating spectroscopic maps with topography, one has to recall that in constant- 
current STM mode a local DOS variation will lead to a change in the tip-sample distance 
z. This can misleadingly induce phantom topographic features. Because the tunnel current 
dependence is exponential on z while only linear on DOS variations, this contribution is 
however modest. We have carefully quantified it to be about 6 pm, that is, more than one 
order of magnitude smaller than the effects described here, see Supp. Info, for details. 

With the above analysis, we arrive to the main result of this Letter: disordered graphene 
on a metallic substrate displays a strong local correlation between local doping and topog- 
raphy. It is important to note that the sign of the local cross correlations x(0) is positive, 
that is, the p-doping character increases with increasing topographic height. Intuitively, 
one expects an infinite graphene sheet above an iridium plane to be p-doped at large dis- 
tances, according to the difference in work functions of both materials. As both surfaces 
come closer, calculations show that the finite electronic wave function overlap leads to a 
cross-over towards ra-type charge transfer [26]. There is thus a wide range of graphene-metal 
distances z at which dEj^jdz > 0, which is here refiected by the sign of x(0)- A contribution 
of curvature effects [5] could also be anticipated. In the present system, the theoretically 
expected contribution of this effect is however nearly two orders of magnitude below the 
observed amplitude of the doping disorder [25] . 

In order to obtain a more microscopic understanding of the role of the local electrostatic 
environment as electron scattering centers, we have tracked the wave vector k distribution of 
interference patterns in the Fourier transforms of DOS maps. The evolution of the dominant 
k, relative to the Brillouin zone corner as a function of the energy level of each map then 
reflects the dispersion relation of the scattered electrons. In graphene, both the large k 
intervalley and the small k intravalley scatterings have been observed by this technique [27] . 
For intravalley scattering, the wave vector transfer q links two points of the circle resulting 
from the intersection of a given Dirac cone and a constant energy plane. A particular 
characteristic of graphene is the suppression of low-energy backscattering [2B] • One therefore 



expects a continuous distribution of q between and 2kF for intravalley scattering, seen as 
a disk in reciprocal space. 

In our data, we observe the DOS maps at energies far from E'^ to display clearly resolved 
structures at length scales smaller than the typical puddle size (Figs. pk-f). Further, the size 
of the observed features decreases with increasing \E — E%\. The maps' Fourier transform 
maps display a disk-like structure (Fig. |3^, see also Supp. Info.) from which we extract 
k = Qmax/'^ at each energy (Fig. jsji, see also Methods). The linear behavior of k{E — E'^) 
evidences the scattering mechanism at work as intravalley. It is actually notable that this 
can be observed in graphene on a metallic substrate at all. 

A fit by the linear dispersion relation of graphene E = E^) ± hvpk (Fig. pp) yields 
vp = 0.89 ± 0.04 X 10^ m/s. In bulk unhybridized graphene, vp is a sensitive probe of the 
strength of electron-electron interactions [21 [29], which can be screened by a large dielectric 
constant e environment. While vp = 2.5 x 10^ m/s on a low-e quartz substrate [SO], it 
decreases to about 1.1 — 1.4 x 10^ m/s on Si02 [3, EI], with a limiting value of 0.85 x 10^ 
m/s expected for infinite screening [30] • The present system can actually be envisaged 
as graphene on a dielectric substrate with divergent e. The experimental value of vp is 
in excellent agreement with the assumption of very effective screening, as expected in the 
immediate vicinity of a metallic ground plane. 

With the above analysis it is seen that (i) the linear dispersion relation of graphene sur- 
vives in the present system at both positive and negative energies and (ii) electron-electron 
interactions in the graphene sheet are strongly screened. This confirms that the present 
system thus behaves in many ways much more like graphene on a dielectric substrate rather 
than a graphene-metal hybridized electronic system, but within a fully screened environ- 
ment. 

One remaining open question concerns the physical origin of the observed quasi-particle 
scattering. In contrast with Si02, coulombic scattering centers in the form of charged impu- 
rities below the graphene do not contribute here, as these should be completely screened out 
by the metallic substrate. It is then tempting to envisage a localization-type interpretation 
to explain the DOS features in real space in Fig. M [31]. We have studied the variations 
of the features studied above (puddles and interference patterns) as a function of magnetic 
field and found no significant change of neither the Ep, nor the G maps up to the maximum 
experimentally accessible magnetic field of 2 T (see Supp. Info.). This leads us to infer that 



the observed interferences stem from long-range scattering by the doping inhomogeneities. 
The present discussion demonstrates the reahzation of strongly screened graphene. The 
two main original features with respect to graphene on a dielectric substrate reside in the 
origin of doping inhomogeneities - set by the graphene-to-metal distance - and the strength 
of electron-electron interactions screening. This is of particular interest for possible quan- 
titative comparisons of the Dirac physics in graphene to quantum electrodynamics (QED) 
[29] . While formally similar, one major difference between graphene and QED is set by the 
relative strengths of the coupling constants: a oc [evp)"^ ^ 1 in graphene and a ~ 1/137 
in QED. The tuning of a by strong screening in a large dielectric constant environment 
would allow recovering the weak-coupling regime in graphene and driving analogies with 
high-energy physics quantitative. 
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Methods 

The experimental set-up is a home made scanning probe microscope operating inside 
a dilution refrigerator. Local tunnel conductance data were obtained using the lock-in 
technique {Vac = 6 mV, /=407 Hz). Prior to cool-down, the sample is heated to 70°C 
overnight while pumping the chamber. During the cool-down, the sample is continuously 
heated well above all other inner parts of the cryostat in order to avoid cryosorption of 
residual gases on the sample. 

The two-dimensional mapping of x(r) yields a sharp peak, at a value fq that is usually 
offset from (0,0) by a few nanometers. This offset is due to drift of the residual lateral relative 
tip-sample position over typical DOS mapping times exceeding 50 hours. The topographic 
maps are therefore shifted by Tq as to maximize x(0)- 



The radius qmax/'^ of the interference patterns in Fig. [3]is defined as the inflection point of 
the angular averaged Fourier intensity (e.g. Fig. Ist), that is, the minimum of its smoothed 
derivative for 0.25 nm^^ < k <1 nm^^. This criterion produces the k{E) dispersion relations 
with least noise. This criterion however overestimates the wave vectors involved by adding 
a constant shift to the detected values of g, which results in the linear dispersion bands 
crossing at A; > in Fig. |3)i. 



FIG. 1: Topography and density of states of graphene on iridium, (a) Scanning tunneling 
microscopy image on graphene on Ir(lll). The image size is 400x400 nm^, tunnel current 1=1 
nA, bias voltage V = 0.57 V). The atomic step in epitaxial Ir(lll) covered by graphene is 2.5 
A high. Inset: zoom-in (2.3x2.3 nm^) atomic resolution image, 1=1 nA, y=0.01V. (b) Local 
tunneling spectroscopy GiV) = dl/dV (see Methods). The Dirac point [eV = Ed, arrow) is 
defined by the minimum of G. 

FIG. 2: Spatial fluctuations of the Dirac point, (a) Dirac point map (color code) superimposed 
with a 3D plot of the long- wavelength topography. Image of 250x250 nm^. (b) Carrier density 
distribution extracted from (a) (see text), (c) Angular average x(|i"l) of the cross-correlation 
function (see text) between the above Ed and topography maps. 

FIG. 3: Spatial inhomogeneities of the local density of states, (a) - (f) G{r,E) maps for 
E = —185; 15; 125; 345; 555 and 700 meV respectively. The color scale is set to cover fluctuations 
of ±25 % with respect to the average G at a given E. As the Dirac point is approached, the 
islands size can no longer be resolved, (g) Power spectral density from angular averaging of Fourier 
transform of (b) (blue) and (c) (green). The inset shows the Fourier transform of (b). The scale 
bar is 0.5 nm^-*^. (h) Energy - wave vector relation extracted from Fourier analysis of the density 
of states maps at energy E. The dashed lines are fits to the linear dispersion relation of graphene, 
yielding vp = 0.89 x 10^ m/s and E% = 0.38 eV. 
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